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Abstract 

The double-differential inclusive di-jet cross section in photoproduction processes 
is measured with the HI detector at HERA. The cross section is determined as 
a function of the average transverse jet energy E^ ts for ranges of the fractional 
energy x J 7 ets of the parton from the photon side. An effective leading order parton 
distribution in the photon is determined at large parton fractional energies for scales 
between 80 < p\ < 1250 GeV 2 . The measurement is compatible with the logarithmic 
scale dependence that is predicted by perturbative QCD. 
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1 Introduction 



In lepton-proton collisions at HERA the cross section is dominated by processes in which the 
lepton radiates a quasi-real photon which then interacts with the proton. In a small fraction 
of such events, jets are formed [0, H), ||, ||, HJ. In the leading-order interpretation of QCD, the 
jets result either from direct photon interactions with the partons of the proton (Fig. |l|a) or 
from resolved photon-proton processes (Fig. |l]b), where the photon interacts with the partons 
of the proton via its own parton content. For resolved photon processes, parton densities of the 
photon /.;/ 7 (x 7 ,/i 2 ) are defined which give the probability of finding a parton i in the photon 
carrying a fraction x 7 of the photon's energy. The scale, at which the photon is probed in the 
parton-parton collision is here considered to be the transverse momentum pt of the scattered 
partons. 

In this paper, a measurement of the inclusive di-jet cross section d 2 a/(dxJ 7 ets d log(£^ ts ) 2 ) 
is presented as a function of the observed parton fractional energy 0.1 < x^ ts < 1 and of the 
squared average jet transverse energy 100 < (Ej! ts ) 2 < 2500 GeV 2 . Throughout the paper we 
denote observed hadronic variables that were corrected for detector effects with the suffix 'jet', 
reconstructed quantities from the detector with the suffix 'rec', and leave variables referring to 
the partons without an additional suffix. All transverse energies are determined with respect to 
the ep beam axis. 

The measurement is compared with QCD calculations of the cross section which require 
parton distributions of the proton and the photon as input. In the kinematic range relevant to 
this analysis, the parton distributions of the proton are well measured. The largest uncertainty 
in the results of the calculations comes from the quark and gluon densities of the photon. The 
quark density functions used in the calculations were determined from measurements of the 
photon structure function in deep inelastic electron-photon scattering at e + e~ colliders. 
These experiments covered the kinematic range 0.05 < x 7 < 0.9 and /r 2 = Q 2 < 400 GeV 2 where 
Q 2 is the virtuality of the photon which probes the quasi- real photon (PETRA |, % |, |H| 0, 
PEP g], TRISTAN 0, [HJ and LEP 0, |l§). In contrast to these F 2 7 measurements, the 
photoproduction of jets is sensitive to both the quark and gluon content of the photon already 
in leading order. A first measurement of the gluon distribution in the photon has been made in 
the range 0.04 < x 7 < 1 at the scale p\ = 75 GeV 2 The comparison of the measured di-jet 
cross section with the QCD calculations provides both a test of the universality of the parton 
densities in the photon and gives new information on the parton densities. 

In addition to predicting the parton-parton scattering processes between the photon and 
the proton, perturbative QCD predicts characteristic features of the quark density / q / 7 in the 
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Figure 1: Diagrams for the di-jet production in ^p scattering: Example for direct (a) and 
resolved (b) photon processes. 



photon [16|. The general characteristics are given by the DGLAP evolution equation for the 
photon 



d/ q/7 (x 7 ,/x) ri 
ln(^/A^ CD ) " aa ^^ j + bas i 7 



d£ 



P 



qq 



/q /7 (z,/i) + P qg 



(1) 



where a and 6 are constant, a and a s are the fine structure constant and the strong coupling 
constant respectively, Aqcd is the QCD scale parameter, and Pij are the splitting functions. 
While the form of the integral on the right hand side is the same as that for hadrons, the 
inhomogeneous term a a P q7 is peculiar to the photon: it reflects the contribution of quarks from 
the pointlike coupling of the photon to a quark-anti-quark pair. Solving (|l]), the asymptotic 
behaviour of / q / 7 at large scale \i and medium x 7 is predicted to be [|l^] 



/q/7 ~ 111 T2 



(2) 



QCD 



The increase of the parton density with increasing scale is clearly different from the behaviour 
of the parton density distributions of hadrons and is referred to as the anomalous hadronic 
structure of the photon. 

In this analysis, the QCD prediction for the scale dependence of the quark densities in 
the photon is tested using the di-jet measurements. Since in the photoproduction of jets the 
incoming quarks and gluons cannot be distinguished, an effective parton density / 7 of the photon 
is extracted from the data. This contains a combination of the quark and gluon densities and 
is valid in leading-order QCD [17|. The measured effective parton density is compared with 
different predictions. 



2 Detector Description 

A detailed description of the HI detector can be found elsewhere |1S[ |. Here we describe only 
those components which are used in the analysis. The HI central tracking system is mounted 
concentrically around the beam-line and covers polar angles between 20° < 8 < 160°. The 
angle 8 is measured with respect to the proton beam direction. The plane perpendicular to the 
z axis is termed r — <j) plane. Measurements of the charge and momenta of charged particles 



are provided by two coaxial cylindrical drift chambers (central jet chambers, CJC) [19|. At 
two radial positions are placed a drift chamber, which provides accurate measurement of the 
z coordinate of charged tracks, and a multiwire proportional chamber (MWPC), which allows 
triggering on those tracks. These are located within the inner CJC and between the two jet 
chambers. The central tracking system is complemented by a forward tracking system which 
covers polar angles 7° ^ 8 ^ 25°. In the present analysis the tracking detectors are used to define 
the vertex position along the beam axis and to support the measurement of the hadronic energy 
flow. 

The tracking system is surrounded by a highly segmented liquid argon (LAr) sampling 
calorimeter [^0j with an inner electromagnetic section consisting of lead absorber plates with 
a total depth of 20 to 30 radiation lengths and an outer hadronic section with steel absorber 
plates. The LAr calorimeter covers polar angles between 4° and 153° and its total depth is 4.5 
to 8 interaction lengths, depending on the polar angle. The region 151° ^ 8 ^ 176° is covered by 
a lead scintillator calorimeter (BEMC). 



A magnetic field of 1.15 T is produced by a superconducting solenoid surrounding the LAr 
calorimeter. The iron flux return yoke surrounding the superconducting solenoid is instrumented 



with limited streamer tubes to provide muon identification and measurement of energy leaking 
from the LAr and BEMC calorimeters. 

All calorimeters provide measurement of the hadronic energy flow. Here, jet measurements 
are restricted to the range which is covered by the LAr calorimeter. 

The luminosity is measured using the radiative process ep — * epj where the photon is 
detected in a luminosity monitor. 



3 Data Selection 

The data were taken in 1994 with the HI detector operating at the lepton-proton collider HERA, 
where positrons of 27.5 GeV collide with protons of 820 GeV. The integrated luminosity used 
for this analysis is 2.9pb _1 . 

The photoproduction events were selected by requiring that the scattered electron remains 
in the beam pipe, undetected in the main part of the HI detector (so-called untagged events). 
The photon virtuality is therefore restricted to Q 2 < 4 GeV 2 , where 80% of the events have 
Q 2 < 0.1 GeV 2 . 

The main trigger for untagged jet events is a summed transverse energy in the LAr calorime- 
ter of at least 6 GeV with the additional requirement of a single energetic cluster of at least 
2 GeV which is validated by a local coincidence of a track originating from the interaction ver- 
tex. Together with triggers that are purely based on the tracking detectors the untagged jet 
events are triggered with at least 80% efficiency over the whole kinematic range. 

The fraction of the lepton beam energy E c that is carried by the photon is denoted by 
y = Ery/E e . It was approximately determined as y rcc from all reconstructed objects in the 
central detector, i.e. calorimeter clusters supported by tracking information, with transverse 
energy an( ^ pseudo-rapidity r^ ec = — In (tan 6/2) using 

y TSC = E Ef c h e~^ c . (3) 
Zhc h 

The reconstructed scaled energy was restricted to the range 0.2 < y Tec < 0.8. The lower cut 
suppresses the contribution from background events that do not result from ep interactions. The 
remaining background is at the level of 1%. The upper cut reduces events from deep inelastic 
scattering at Q 2 > 4 GeV 2 , where the scattered electron remained unidentified, to less than 1% 
of the total sample. 

The jets were found using a cone algorithm [^] with cone size R = 0.7 and jet pseudo- 
rapidities in the range — 0.5 < rf ec < 2.5 in the HERA laboratory system. In this implementation 
of the cone algorithm, the jet cone position is chosen such that the transverse energy of the jet 
is maximal. The jet with the highest transverse energy is constructed first, the jet finding is 
then continued using only the remaining energy flow outside of the first jet cone. In case of 
overlapping jet cones, all energy flow in the overlap region is assigned to the jet with the higher 
transverse energy. 

After these selection cuts, the total number of events with at least two jets reconstructed 
with a transverse energy above 7 GeV was 6499. 



4 QCD Calculation of the Di-Jet Cross Section 



In leading-order QCD the differential ep cross section for di-parton production can be written 
as follows: 

d4 ^ _ 1 fjy/M y /i/ 7 QE 7 ,Pr) /j/pOWr) i M ..( cos fl*)|2 _ (4) 

dy dx 7 dx p dcos 9* 327rs ep y ^rr 1 x 1 x p %3 

In this expression, s ep denotes the squared lepton-proton center-of-mass energy available at 
HERA, ^/s cp = 300 GeV. The flux of photons radiated off the electron with fractional energy y 
is predicted by QED and is denoted by / 7 / e . The fractional energy of the parton in the photon is 
given by rr 7 and the parton density function of parton i in the photon by The corresponding 
variables for the proton are x p and fj/ p - Note that, in contrast to the quark distribution in the 
photon mentioned in (jl]), here represents all components in the photon including the direct 
photon processes and the quarks and gluons for resolved photon processes. The last term in (||) 
contains the distribution of the parton scattering angle 9* in the parton-parton center-of-mass 
system in the form of the QCD matrix elements Mij(cos9*). 

The parton transverse momentum px, which is here identified with the scale of the process, 
can be expressed in terms of the four observables in (Q): 

Pt = J s ep V x 1 x p sin 2 9* . (5) 
For the correction of the data, and for comparison to the measured jet cross sections, the 



two event generators PYTHIA [£2] and PHOJET [26, 27 were used. In both generators jet 



production is based on the leading-order di-parton cross section as given in (0). 



The PYTHIA 5.7 event generator [22| was used to simulate photon-proton interactions. As 
well as the leading-order cross section (Q), PYTHIA includes initial- and final-state parton radi- 
ation effects which are calculated in the leading logarithmic approximation. The strong coupling 
constant a s is calculated in first order QCD using Aqcd = 200 MeV for 4 quark flavours. Mul- 
tiple parton interactions were generated in addition to the primary parton-parton scattering. 
They are calculated as leading-order QCD processes between partons from the photon and pro- 
ton remnants. The transverse momentum of all parton-parton interactions was required to be 
above px = 1.2 GeV. This cut-off value has been found to give an optimal description of the 
transverse energy flow outside the jets Q . A good description of this energy flow is important 
since the contributions of multiple parton interactions inside the jet cones alter the jet rates 
considerably. For tests of the uncertainties in the determination of the effective parton distribu- 
tion of the photon arising from the contribution of multiple parton interactions, an event sample 
was generated with px > 1.4 GeV, which gives an energy flow that is significantly smaller than 
that observed in the data. To assess the significance of higher-order effects, another sample was 
generated without initial-state parton showers. GRV-LO parameterizations were used through- 
out for the parton distributions of the photon and the proton p3| , |24| a nd hadronization was 
modelled with the LUND string fragmentation scheme (JETSET 7.4 |25fl). 



The PHOJET 1.06 event generator p(l, 27] is based on the two-component Dual Parton 



Model |2q| . It calculates parton-parton scattering using (f|). Again, Aqcd is 200 MeV for 4 
quark flavours. PHOJET incorporates very detailed simulations of both multiple soft and hard 
parton interactions on the basis of a unitarization scheme . It includes initial- and final-state 
hard parton radiation effects. The lower momentum cut-off for hard parton interactions was set 
to pt = 2.5 GeV. In contrast to the PYTHIA generator, small variations of this cut-off value do 
not have a large influence on the results obtained by PHOJET, due to the unitarization scheme. 
For the fragmentation the LUND string concept is applied (JETSET 7.4 pi). The GRV-LO 



parton distribution functions for the photon and the proton were also used here. The effects of 
higher-order corrections were checked by generating a sample of events without hard initial-state 
parton showers. 

For comparisons with the measured jet cross sections, analytic QCD calculations were pro- 
vided by ref. |3(J. These calculations include the next-to-leading order (NLO) QCD matrix 
elements. They use the NLO parton distribution functions GRV-HO |23| and GS96 [31] for the 
photon and the CTEQ-4M distributions [32] for the partons of the proton. 



5 Measurement of the Di-Jet Cross Section 

In order to study observables closely related to the parton distributions in the photon, the 
differential di-jet cross sections were analysed in terms of the variables E^ ts and x^ ts . The 
observable E^ ts was calculated using the two jets with the highest E^ in the event 

pjetl . p jet2 

4 ets = (6) 

and was required to be above E^ ts = 10 GeV. 

The ratio of the difference and the sum of the transverse energies of the jets was required to 
be in the range 

I F jetl _ T-,jet2, 

' " i^<0-25. (7) 



+ E% t2 



The cuts (|) and (@) ensure that the transverse energy of all jets is above E^ = 7.5 GeV, without 
using the same E^ cut-off for both jets. In next-to-leading order analytical calculations of the 
di-jet cross section this avoids introducing a dependence of the result on unphysical parameters. 

The observable x^ ts was calculated using the transverse energies and pseudo-rapidities of 
the two jets and the photon energy: 

T jets _ T e 
7 " 2yE c ■ W 

Here, the interval 0.1 < x^ ts < 1 is considered. 

The average rapidity of the two jets was required to be in the range 

rjjetl I „jet2 

< £ P < 2 . (9) 

This variable is approximately the rapidity of the parton-parton center-of-mass system in the 
laboratory frame of reference. The difference in the jet pseudo-rapidities was required to be 
within 

\Ari cts \ < 1 (10) 

which corresponds to | cos#*'j ets | < 0.46. Together with @ this cut ensures that all jets are in 
a region in which the hadronic energy is well measured in the detector. 

In Fig. [2] the ep double-differential di-jet cross section is shown as a function of the transverse 
energy Ej, s in bins of the fractional parton energy x^ ts . The event kinematics corresponds to 
corrected photon fractional energies in the range 



0.2 < y < 0.83 



(11) 
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Figure 2: The measured inclusive di-jet ep cross section is shown as a function of the squared 
jet transverse energy E^ ts for ranges of the reconstructed parton fractional energy x^ ts (full 
circles). E } ^ ts is the average transverse energy of the two jets with the highest in the 
event. The average pseudo-rapidities of the jets were between < (r^ etl + r^ * 2 ) /2 < 2 and their 
difference was below IAt^ ** 3 ! = 1. The transverse energies of the jets were restricted to the range 
\E^ tl — E^ t2 1 / (Eji* + E^ t2 ) < 0.25. The cross section is integrated over the photon virtuality 
Q 2 < 4 GeV 2 and the relative photon energy 0.2 < y < 0.83. The inner error bars represent 
the statistical errors, the outer error bars give the statistical and systematic errors, added in 
quadrature. The data are compared to the QCD simulation of the PYTHIA generator using the 
GRV-LO parton distribution functions (dashed curve) and to analytical next-to-leading (NLO) 
QCD calculations of the di-jet cross section |?0| using the GRV-HO (full line) and GS96 (dotted 
line) photon parton distributions. 
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Table 1: Inclusive double-differential di-jet cross section in bins of x J 7 ets and log(.Ep ts /GeV) 2 . 
The limits of the bins are listed in the first four columns. For each bin a cross section number 
is given together with the statistical and systematic uncertainties. The error induced by the 4% 
uncertainty of the hadronic calibration of the calorimeter (E-scale) is listed separately. 



and photon virtualities of 

Q 2 <4GeV 2 . (12) 



The di-jet cross section was corrected for detector effects using an unfolding procedure |33[. The 
procedure requires a 4-dimensional correlation function between the variables reconstructed 
in the detector (x™ c ,Ej? c ) and the corrected observables referring to the hadronic final state 
(x^ ts , E^ ts ). This correlation was calculated using the PYTHIA generator together with a de- 
tailed simulation of the detector. The procedure modifies the di-jet cross section of the PYTHIA 
calculation without detector effects so as to ensure that the predicted di-jet distributions as a 
function of (x 7 ec , E^ c ) agree with the data. Further details of the correction procedure can be 
found in j34|. The inner error bars reflect the statistical errors. The outer error bars shown 
are the statistical and systematic errors added in quadrature. The main contribution to the 
systematic error results from a 4% uncertainty in the knowledge of the calorimeter hadronic 
energy scale (Table [l]) . Variation of the unfolding parameters affects the measured cross section 
by less than 10%. The uncertainty in the determination of the trigger efficiency using redundant 
triggers is below 10%. The luminosity measurement is accurate to 1.5%. 

The measured cross section is compared to the simulation of the PYTHIA generator (dashed 
curve) using the GRV-LO parton distribution functions for the photon and the proton p3| |24| ], 
This calculation includes the leading-order QCD matrix elements together with higher-order 
QCD effects which are simulated by a parton shower mechanism, multiple parton interactions 
and hadronization. All x]^ ranges are dominated by resolved photon processes, except for the 
highest region, x J 7 cts > 0.75, where the direct photon processes give the dominant contribution. 
The PHOJET generator, using the same parton distribution functions and the same choice of 
the scale, gives results that are compatible with the PYTHIA calculation to within 20% (not 
shown). Also shown are two analytical calculations of the di-jet cross section in next-to- leading 
order QCD @ which use the GRV-HO || (full curve) and GS96 || (dotted curve) photon 
parton distribution functions. In both calculations the CTEQ-4M |3^] parameterization of the 
proton parton distribution is used. In these calculations the jet finding is based on the three 
final state partons. Since fragmentation effects and underlying event energy from multiple parton 
interactions are not included, deviations from the measured cross sections are expected. Multiple 
parton interactions increase the jet energy and therefore the jet cross section. The influence on 
the cross section has been estimated to be below 40% at small x^ ts ~ 0.1 and simultaneously 
small M? ts , and is negligible at xl^ ts ~ 1 or E 3 ^ ts > 15GeV. The influence of the fragmentation 



r jets 



process is maximal in the region of low Ej, and raises the jet cross section by at most 20% [35]. 



All calculations give a satisfactory description of the data overall, except for some deficiencies 
in the lowest x^ ts range and the two highest x^ ts ranges. In these regions, the resolved photon 
contribution is not well constrained from measurements of the photon structure function F% at 
e + e~ colliders. The two GRV-LO and the GS96 parton distribution functions each give a good 
description of the F^ measurements. The successful application of these parton distributions in 
the calculation of the di-jet cross section in ep scattering supports the universality of the photon 
parton distributions. However, in detail, the differences between the measurement and these 
calculations show that the jet data give new constraints on the parton distribution functions, 
especially in the region at large x^ ts and high E^ ts . 



6 Effective Parton Distribution of the Photon 



In order to extract the parton distributions of the photon from the data, /j/ 7 in (Q) needs to 
be disentangled from the sum over the different parton processes. This is not possible without 
further assumptions. 



Here we follow the approximation procedure developed in [17]. The method is based on 
the leading-order matrix elements for the partonic scattering processes qq 1 — * qq', qg — * qg, 
and gg — * gg which give the dominant contribution to the total resolved photon di-jet cross 
section. The shapes of the angular distributions of the scattered partons are similar for all of 
these processes. The squared matrix elements differ in rate by the ratio of the colour factors 
Ca/Cf = 9/4. The approximation uses only one Single Effective Subprocess matrix element 
Mses an d combines the quark and gluon densities into new effective parton density functions 

/ \ 9 

/ 7 (x 7 ,p T ) = (/qAyOWT) +/q/ 7 ( X 7'PT)J + jfgh(x-y,PT) and ( 13 ) 

iif 

/ \ 9 

fp(Xp,P T ) = Y (/q/pOWT) + h/p(Xp,PT)) + T /g/p(«p,Pr) • ( 14 ) 

The sums run over the quark flavours. The product / 7 / p |Mses| 2 then replaces the resolved 
photon contributions to the sum ^ - fi^fjj-^Mi^ 2 of equation (||): 

^ ~ ^ AM) fel) |M SES (cosr)| 2 . (15) 



dy dx 7 dx p dcos 6* 327rs ep y x y x p 

Here / 7 (x 7 ,p|i) can be directly determined from the measurement of the di-parton cross section. 
In the kinematic range which is relevant to this analysis the accuracy of the 9/4 weight in (|l3| ) 
was checked with a PYTHIA calculation of the di-parton cross section based on (^|) . The relative 
contributions of the quark and gluon induced di-jet events agree with this weight to better than 
5% in the range 0.2 < x 7 < 0.7 considered here. 



6.1 Test of the Factorization of the Di-Jet Cross Section 

Before the effective parton density / 7 of the photon is extracted, it is shown that the factorization 
of the di-jet cross section ( |l5| ) into effective parton densities is compatible with the measured 
data. If the factorization approximation is good, the shape of the parton fractional energy dis- 
tribution xp ec from the proton side should be independent of the parton fractional energy from 
the photon side x™ c . In order to test this hypothesis, the four dimensions of the reconstructed 
di-jet distribution N(y* ec , x 7 ec , x p ec , #* ,rcc ) are reduced to two dimensions, using the two inde- 
pendent observables x p ec and y rec by applying the requirements #*> rcc 90° (|Ar/ rcc | < 1) and 
yTec x rcc ^ constant (0.1 < y rcc x™ c < 0.2). The latter cut implies that the energy which enters 
the parton-parton scattering process from the electron side is constant. x p ec is calculated from 
the reconstructed jet energies and pseudo-rapidities using Xp 6C = e^° c + E^ 2 e^") / (2E P ). 

In Fig. ||, the x p ec distribution is shown in 4 different y iec intervals. All distributions were nor- 
malized to the data sample covering the full y mc interval, 0.2 < y vcc < 0.8. The measured ratios 
are compatible with having the same shape in all four y Tec bins. 

Since the average parton fractional energy {x™ c ) of the four x p ec distributions varies between 
(x 7 ec ) = 0.22 and (x 7 ec ) = 0.69, the observed x T p ec distributions are also independent of x™ c . 
Therefore, within the precision of the data, factorization holds in (pl|) and so the ansatz used 
to extract an effective parton distribution for the photon is meaningful. 



6.2 Effects Influencing the Correction to Leading Order Parton Variables 

In order to correct from the jet observables to the parton variables three effects have to be taken 
into account. These are multiple parton interactions, higher-order QCD effects and fragmenta- 
tion effects. 
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Figure 3: The distributions of the reconstructed fractional energy x p ec of the parton from the 
proton side are shown in four bins of the photon fractional energy y rec . All distributions were 
normalized to the x™ c distribution of all events in the whole interval 0.2 < y TCC < 0.8, diVo/dXp CC . 
Only statistical errors are shown. The dotted line represents a fit to the ratio assuming a flat 
distribution. 



In the generators used here, the fragmentation is simulated using the well tuned Lund string 
fragmentation model (JETSET 7.4). 

Multiple parton interactions lead to an additional energy flow in the event which affects the 
jet rates. These effects were studied in a previous analysis using the transverse energy flow 
outside the jets and energy-energy correlations. These are well simulated, for example, by the 
PYTHIA generator g]. 

Higher-order QCD effects can be studied by looking at multi-jet production. The relative 
contribution of events with three or more reconstructed jets above a transverse jet energy of 
7GeV amounts to 8% of the total event sample. This contribution is well described by the 
PYTHIA simulation. 

Further results of higher-order QCD effects are an imbalance between the transverse energies 
of the two highest transverse energy jets and a deviation from a back-to-back configuration in the 
azimuth. For this study, the region of reconstructed parton fractional energies above x^ ec = 0.4 
is considered in order to minimize contributions of multiple parton interaction effects []|]. The 
missing total transverse energy in photoproduction events is small and was here required to be 
below -E^^iss = ^ GeV to ensure that the transverse jet energies were well measured. 



In Fig. |]a, the shape of the transverse energy imbalance AE^p c = \E^\ — E^\ between the 
two jets with the highest E^ c is shown. Fig. |^b shows the azimuthal difference between the two 
jets A(J) rec = \4> r i c — </>2 ec |- The shape of the AE^ C and Ac/> rec distributions is described by a 
PYTHIA calculation which includes hard initial-state parton showers (full histogram in Fig. 
This calculation also gives a good description of the AE™ c and A(p rec distributions at x™ c < 0.4 
(not shown). The dashed histogram represents a calculation of the PHOJET generator in a 
version that does not include hard initial-state parton radiation effects: this calculation gives 
too small AE^p c and too large A(j) Tec . We conclude that higher-order QCD effects are well 
modelled by parton showers. 
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Figure 4: a) The shape of the distribution of the uncorrected jet transverse energy imbalance 
AE?p c = \ETp\ — -Ep^l i s shown. The distribution is normalized to the maximum number of 
entries iV max in a single bin. The reconstructed parton fractional energy was required to be 
x™ c > 0.4. The full histogram shows the prediction of the PYTHIA generator which includes 
hard initial-state parton radiation effects (ISR). The dashed histogram shows the calculation of 
the PHOJET generator in a version without hard initial-state parton radiation (no ISR). Both 
calculations include a detailed simulation of the detector effects, b) The shape of the uncorrected 
distribution of the azimuthal difference between the two jets Ac/> rec = |^ ec — (/>?> ec | is shown. The 
histogram assignments are as in a). 



6.3 Extraction of the Effective Parton Distribution Function 

For the extraction of the effective photon parton distribution function, / 7 , the measured di-jet 
cross section was corrected to the level of the leading-order di-parton cross section using the 



same unfolding procedure |33[ as applied in the analysis of the jet cross section described above. 
In this case the correlations between (x^ ts ,E^ ts ) and the parton variables (x 7 ,pr) are used to 
correct for the fragmentation, higher-order QCD effects and underlying event energy effects. 
As discussed in the previous paragraph, all these effects are well modelled by the PYTHIA 
generator. The model dependence of the corrections was checked using the different generated 
event samples described in Sect. [|. Comparison of the di-parton cross section from data and a 
PYTHIA calculation using (||) together with the GRV-LO parton distributions and with px as 
the scale then gives the effective parton distribution of the photon in the data: 

~ data _ -GRV-LO d 2 (j DATA /(da; 7 d \ogPr) 

jl ~ /7 X d 2 cT PYTHIA ' GRV - LO /(da; 7 d log p 2 T ) ' U 



In Fig. [|the measured effective parton distribution of the photon and its scale dependence are 
presented. The data are shown in two intervals of the parton fractional energy a) 0.2 < x^ < 0.4 
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Table 2: The effective photon parton distribution x 7 / 7 is shown as a function of the parton 
fractional energy x 7 and the transverse momentum p\ of the scattered partons. The systematic 
errors are separated into contributions from experimental uncertainties and model dependencies. 
The total error (<7 to tai) corresponds to the quadratic sum of the statistical (<7 s tat) and systematic 
(o- S y St ) errors. 

and b) 0.4 < x 7 < 0.7. The inner error bars represent the statistical errors, the outer error 
bars the quadratic sum of the statistical and systematic errors. In addition to the errors of the 
di-jet cross section, which were described in Sect. |5| the error bars include the uncertainty in 
the correction resulting from multiple parton scattering effects which was determined using the 
two event generators PYTHIA and PHOJET (Table |). 

For comparisons with the data, the effective parton distributions of the photon and the pion 
were calculated from the GRV-LO parameterizations p4| , [36fl . The pion parton distribution 
was scaled by the factor k ■ Aira/ f 2 ~ 0.9 a where /p/47r ~ 2.2 represents the probability that a 
photon converts into a p meson and k ~ 2 effectively accounts for contributions of heavier vector 
mesons |23| . This parton distribution therefore gives the vector meson dominance picture of the 
photon and is shown as the dashed curve (VDM) in Fig. 0. It differs from the measurement 
both in shape and absolute rate. 

The GRV-LO parameterization of the photon parton distribution functions was chosen to 
resemble that of a hadron at small scales, [i ~ 0.5 GeV (not shown). At large scales /i>l GeV, 
the parton densities of the photon reproduce the measurements of the photon structure function 
which is observed to rise with fi 2 = Q 2 at a rate which is compatible with the logarithmic 
increase that results from the pointlike term in the DGLAP QCD evolution equation (|l]) (see 
e.g. @). 

The measurements made here show, for the first time in the context of photoproduction, 
the above-mentioned logarithmic dependence of photon structure on the scale, pt, at which 
that structure is probed. The GRV-LO parameterization of photon structure (full curve in 
Fig. |5|), which includes the effects of the pointlike term in the DGLAP equations, describes the 
measurement. If this term is excluded, as in the VDM model calculations shown here (dashed 
curve), it is not possible to describe the data. Further, as information on the photon quark 
distributions is obtained from the i 7 ^ measurements, it is possible to identify the contribution 
they make to the effective parton distribution measured here (dotted curve). The difference 
between this and the measured values is a consequence of the gluon content of the photon. 
This is observed to contribute about 20% to the total effective parton distribution in the range 
0.4 < x*y < 0.7 and about 50% in the lower 0.2 < x 7 < 0.4. 

The precision of the jet measurement is similar to that obtained in the e + e~ experiments. 
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Figure 5: The leading order effective parton distribution of the photon x 1 f 1 (x 1 ,p\) = 
x 1 (f q / 1 (x 1 ,p^) + 9/4/ g / 7 (x 7 ,p^)) is shown as a function of the squared parton transverse mo- 
mentum p\. The / q / 7 represents here the sum over all quarks and anti-quarks and / g / 7 is the 
gluon density. The distribution was divided by the fine structure constant a and averaged over 
parton fractional energies in the ranges a) 0.2 < x 7 < 0.4 and b) 0.4 < x 7 < 0.7. The inner 
error bars represent the statistical errors, the outer error bars give the statistical and system- 
atic errors, added in quadrature. The data are compared to the effective parton distribution of 
the photon which includes the pointlike coupling of the photon to quarks (full curve) and to a 
vector meson dominance (VDM) ansatz for the photon (dashed curve) . The dotted curve shows 
the quark part of the effective parton distribution. All three curves were calculated using the 
GRV-LO parton distribution functions. 



The data extend the kinematic region where the parton distributions of the photon are measured 
to the scale p 2 T = 1250 GeV 2 . 

7 Summary 

A measurement of the double-differential inclusive di-jet cross section in terms of the parton 
fractional energy and the transverse energy scale was presented from HI photoproduction data. 
This measurement constrains the quark and gluon distributions of the photon with a precision 
that is competitive with two photon data. In addition, a new kinematic range up to x J 7 ets « 1 

and (£^ cts ) 2 = 2500 GeV 2 is covered. 

For the first time, an effective parton distribution of the photon was extracted from the 
data. The observed scale dependence shows an increase with -p\ that is compatible with the 
logarithmic increase predicted by perturbative QCD calculations of the parton content of the 
photon. 
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